Structural phase transitions of vortex matter in an optical lattice 
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We consider the vortex structure of a rapidly rotating trapped atomic Bose-Einstein condensate in 
the presence of a co-rotating periodic optical lattice potential. We observe a rich variety of structural 
phases which reflect the interplay of the vortex-vortex and vortex-lattice interactions. The lattice 
structure is very sensitive to the ratio of vortices to pinning sites and we observe structural phase 
transitions and domain formation as this ratio is varied. 
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Structural phase transitions (SPTs) occur in a wide 
range of physical and chemical systems 0. Studies of 
SPTs are both of technological importance and of fun- 
damental interest. Technologically, an understanding of 
SPTs is essential to the study of many materials; while 
fundamentally SPTs are related to a variety of questions 
in statistical mechanics, crystallography, magnetism, sur- 
face science, to name a few. Among the most intensively 
studied systems related to SPTs are graphite intercala- 
tion compounds, niobates, adsorbed molecular monolay- 
ers and vortex matter in type-// superconductors. In this 
Letter we introduce a new system for SPT study: that 
of vortex matter created in a trapped Bose-Einstein con- 
densate (BEC) formed in the prescence of a co-rotating 
optical lattice (OL). This system has the advantage of be- 
ing experimentally realizable, tunable over a wide range 
of interaction parameters, and describable by an analytic 
ab initio theory. 

We take inspiration from two recent developments in 
atomic BECs. The first is the observation of vortex lat- 
tice excitations in rapidly rotating BECs Q ■ The second 
is the set of proposals and demonstrations concerning pe- 
riodic optical lattice potentials created using holographic 
phase plates or amplitude masks 3 . This latter tech- 
nique means that a rotating OL (at desired angular fre- 
quency) can be realized by simply rotating the OL phase 
plates or masks. The amenability of BECs to imaging, 
experimental control and theoretical description makes 
the atomic condensate an attractive test beds for many 
phenomena frequently encountered, but often difficult to 
study in other systems. One such example is given by 
the recent spectacular realization of the superfluid-Mott 
insulator quantum phase transition in a condensate con- 
fined by an optical lattice [j| following a proposal of 
Jaksch et al. This transition was theoretically pre- 
dicted over a decade ago Q in the context of liquid 4 He 
absorbed in porous media. 

In our system, the phase transition is observed as a 
shift in the periodicity and symmetry of the vortex lat- 
tice as optical potential parameters are tuned. Funda- 
mentally, this situation is analogous to that of a type-// 
superconductor subject to external magnetic field and 
artificial periodic pinning 0, 0, 0] . In the analogy, the 
angular rotation Vtz of the BEC plays the role of the 



magnetic field and the peaks of the OL play the role of 
the pinning centers. In contrast to the pinning lattice 
in typical superconducting samples, the periodicity and 
depth of the optical lattice may be dynamically tuned in 
an atomic system. Furthermore, the clean microscopic 
physics of atomic condensates makes a first principles cal- 
culation possible based on a mean-field treatment. This 
is in contrast to the superconductor case, where theo- 
retical calculations frequently rely on elegant, but phe- 
nomenological models or molecular dynamics simulations 

a 

We work in a pancake shaped geometry. For the BEC, 
this is justified because at high rotation rates, the cen- 
trifugal forces reduce the radial trapping frequency u>± 
and the condensate may be accurately assumed to be 
frozen into the harmonic oscillator ground state in axial 
direction (z). We therefore integrate the axial degree of 
freedom, obtaining an effective two-dimensional system 
with a renormalized coupling constant. In a frame rotat- 
ing with angular velocity Qz the transverse wave function 
4>(x, y) obeys the two-dimensional time-dependent Gross- 
Pitaevskii equation: 



dt 



+ V + Mat + u\^\ 2 - M - m. 



4>. (1) 



Unless otherwise specified we choose ou r units for time, 
length and energy as 1/ oj± , ^H/(muj±) and fruj± , respec- 
tively. Here r 2 = x 2 +y 2 , U is the (renormalized) effective 
2D nonlinear interaction coefficient, \x the chemical po- 
tential, L z the z-component of the angular momentum 
operator, and Vj at = Vb[sin 2 (fcx) + sin 2 (fey)] is the optical 
lattice potential. In this work, we only consider the case 
where the optical lattice is co-rotating with the conden- 
sate, thus Vi a t in Eq. Q is time- independent. 

For each chosen value of the OL periodicity (denoted 
by a = n/k), we vary the peak-to-peak potential ampli- 
tude of the optical lattice Vq from to Vm ax , a value at 
which the lattice is fully pinned. We obtain the ground 
state of the system for a given set of parameters (Vo, a, f2) 
by propagating Eq. imaginary time (i.e. by steepest 
descent). To ensure that the obtained structure corre- 
sponds to the ground state, we usually start with several 
different initial trial wave functions possessing different 
symmetries for each set of parameters. In our calculation, 
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FIG. 1: Structure of the ground state. Each intersection of the dashed lines represents an antinode of the optical potential. 
Filled (•) and empty (o) circles represent the positions of singly and doubly quantized vortices, respectively. Only the positions 
of vortices within the Thomas-Fermi radius are shown. The Thomas- Fermi boundary Rtf = V^P is represented by the circle in 
the upper-left plot. The upper, middle and lower rows correspond to a — 4e/3, e and 2e/3, or 77 «2, 1.15 and 1/2, respectively. 
From left to right, Vo = 0.05, 0.2, 0.5, 2.0 and 5.0, respectively. 



we use the 87 Rb parameters of the JILA experiment 0: 
uj_ = 2tt x 8.3Hz, uj z = 2tt x 5.3Hz, fl = 0.95w_l, and 
the total particle number is taken to be 2 x 10 5 . Evi- 
dently, both the atomic and optical parameters will af- 
fect the ground state. For simplicity, in this work we 
focus on the effect resulting from variation of the OL. 
Our variables are therefore the depth (Vo) and periodic- 
ity (a) of the optical potential. Specifically, a defines the 
density of pinning sites which occur at the peaks of the 
OL potential whereas Q determines the density of vor- 
tices. As found in previous studies |9j, we observe that 
the structure of the fully pinned vortex lattice (FPVL) 
is very sensitive to the ratio i] defined as r] = N v /N p , 
where N v and N p are the density of vortices and pinning 
sites, respectively. For the OL potential we use here, 
N p = l/a 2 . In the absence of the OL, N v = 2/(y / 3e 2 ), 

where e = \J 2ir/ (\/3~fi) is the inter-vortex spacing in a 
completely unpinned (Vo = 0) vortex lattice |ll|. This 
value of N v is not significantly changed by the presence 
of the OL. Hence we take r\ » 2a /(\/3e 2 ). In brief, we 
find that when r) w 1/2, the FPVL has a checkerboard 
structure and when r\ « 1, we obtain a square lattice 
while when r\ w 2 we obtain a square lattice of doubly 



charged vortices. 

The trends in our results are graphically illustrated in 
Fig. n our main result. For three values of a the ground 
state vortex structure is plotted for increasing values of 
Vq. The upper, middle and lower rows correspond to dif- 
ferent pinning site densities represented by a = 4e/3, e 
and 2e/3, respectively. In this situation, the angular mo- 
mentum of the condensate is carried in singly quantized 
vortices which organize into a triangular Abrikosov lat- 
tice due to the logarithmically repulsive inter-vortex in- 
teraction. Such vortex lattices have been observed [l^ . 
perturbed @] and accurately described |l| by several 
groups in recent years. 

As is not surprising, for sufficiently small Vo the vortex 
lattice maintains a triangular geometry with only slight 
distortions (see first column of Fig. In the opposite 
extreme, at sufficiently large values of Vo, all the vor- 
tices are pinned to the antinodes of the OL potential, 
mirroring the geometry set by the OL potential. What 
is remarkable are the states that exist between these ex- 
tremes. At a = 4e/3 (upper row of Fig. QJ, there are 
about twice as many vortices as pinning sites (77 ~ 2), 
and the FPVL is found to be a square lattice of doubly 
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defined as: 
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FIG. 2: Structure factors as functions of Vb for a = 4e/3, e 
and 2e/3. St,s,c = |5(kT,s,c)|, where k T ,s,c = (2n/e)x — 
[27r/(v3e)]y, 2fe?/ and fex — fey, representing respectively one 
of the two fundamental reciprocal vectors for the triangular, 
square and checkerboard lattices. 



quantized vortices which is commensurate with the OL, 
a situation analogous to the second matching field case 
for the superconductor system in which the observation 
of a square double-quanta vortex lattice was reported in 
Ref. 10]. At a higher pinning site density defined by 
a = e (middle row of Fig.^ 77 « 1.15), most of the antin- 
odes of the OL are occupied by singly quantized vortices 
with the exception of three doubly quantized vortices. 
And finally at a = 2e/3 (lower row of Fig. 77 1/2), 
every next-nearest-neighbor site of the OL is occupied by 
a singly quantized vortex and the FPVL forms what can 
be described as a "checkerboard" square lattice rotated 
45° with respect to the OL, once again analogous to the 
half matching field case of the superconductor system . 

For intermediate values of Vq, the vortices form struc- 
tures in between the triangular lattice and the FPVL, 
the details of which also depend on the period of the 
OL. For example, at a = 4e/3, for which the FPVL is 
a square lattice of doubly charged vortices, we observe 
bound pairs centered around the OL pinning sites for 
the potential depth in the range 0.5 < Vq < 4.0. We 
point out that the orientation of each pair is orthogo- 
nal to all adjacent ones. As Vq increases, these pairs are 
more and more tightly bound and eventually all pairs 
collapse onto the corresponding pinning sites thus form- 
ing doubly quantized vortices. Further increasing the OL 
period, vortices with higher and higher winding number 
will start to appear in the FPVL. 

In order to characterize the SPT more quantitatively, 
we calculate the structure factor [9j of the vortex lattice 



where i labels individual vortices, rij and r$ are the 
winding number and position of the i-th vortex, while 
N c = JV rii is the total winding number. For a familiar 
crystal lattice, |5(k)| displays peaks at the correspond- 
ing reciprocal lattice vectors. Here, we focus on follow- 
ing three cases: (1) the triangular Abrikosov lattice in 
the absence of the OL (St), (2) the square (Ss) and (3) 
the checkerboard (5c ) lattices defined by the OL. Each 
lattice structure has two fundamental reciprocal vectors. 
In this instance it is sufficient to calculate the structure 
factor along one dimension, i.e., for one of the reciprocal 
vectors. Our results are displayed in Fig. |2 We observe 
that as Vq is increased from zero, the triangular lattice is 
destroyed over a very small range of Vq as St exhibits a 
sudden jump, indicative of a first order transition which 
is physically expressed by the motion of vortices towards 
the pinning sites. When there are more vortices than 
pinning sites (77 > 1), the surplus vortices get pinned at 
a comparatively "slow" pace as a consequence of repul- 
sion experienced from vortices which are already pinned. 
Conversely, when the number of pinning sites exceeds 
that of vortices (77 < 1), the FPVL is quickly established. 

We now turn to explore the structure of the FPVL 
for various values of 77. In the case 77 is not close to an 
integer or inverse integer value, we observe the coexis- 
tence of sub-lattices of different geometry bounded by 
domain walls which always lie along the diagonal of the 
OL. For -q = 1.8 we have the coexistence of two sub- 
lattices formed by doubly and singly quantized vortices, 
respectively (see Fig. |3Ji). For this FPVL the domains 
have striped parallel walls (represented by the solid lines 
in the figure) , indicating inter-domain repulsion [Til Il5| . 
This can be understood as a consequence the repulsive 
vortex-vortex interaction. For the slightly smaller value 
of 77 = 1.6, these two sub-lattices form two interlock- 
ing checkerboard structure (see Fig. Et>) ■ As a decreases 
further to 0.95e (Fig. |3t), all doubly quantized vortices 
have disappeared and the vortex lattice becomes com- 
pletely commensurate with the OL, with each pinning 
site hosting a singly quantized vortex. This represents 
the fully matching case of 77 = 1. Upon further decrease 
of a or 77, Fig.|21i and e show that more pinning sites be- 
come unoccupied and the checkerboard domain begins to 
cover the whole lattice. In this ground state, the walls be- 
tween the square and checkerboard domains are crossing 
each other, signifying an attractive domain wall inter- 
action which may be intuitively understood as resulting 
from the tendency of the vortices to occupy the vacant 
pinning sites. 

In the study of domain wall formation, we see again 
distinct advantages of BEC vortex lattices over alterna- 
tive systems. Although similar domain formation was ob- 
served earlier in superconducting systems |10L llOj , their 
origin could not be clearly established because of addi- 
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FIG. 3: Structure of the FPVL. From left to right, 77 =1.8, 1.6, 1.0, 0.8 and 0.6, corresponding to a/e =1.25, 1.176, 0.95, 0.833 
and 0.741, respectively. Solid lines represent domain walls. 



tional defects likely to be present in the experimental 
samples. By contrast, in the atomic BEC system, the 
OL provides us with a defect-free periodic pinning poten- 
tial and thus frees us to confidently investigate the more 
fundamental factors controlling the dynamics of domain 
formation and hence the structure of the domain walls. 
Recent success at directly imaging and calculating ^| 
vortex lattice excitations in BECs presents very exciting 
possibilities for studying, in unprecedented detail, the 
dynamics of the FPVLs obtained in this work and opens 
the door to more ambitious studies of SPTs in unconven- 
tional geometries. 

In summary, we have theoretically investigated the vor- 
tex state of a rapidly rotating condensate in a periodic 
optical potential. We have found that the vortex lattice 
exhibits a rich variety of structures depending on the pa- 
rameters of the optical potential. In the future it will 
be interesting to investigate the detailed dynamics (e.g. 
time evolution) of the phase transitions between various 



structures and to relate the transitions to the properties 
of the condensate such as its collective excitation modes. 
Such studies will certainly shed new light on many other 
systems displaying structural phase transitions. In the 
current work, we have assumed that the optical poten- 
tial is co-rotating with the condensate. We also plan to 
generalize this work to the case when the OL and the 
condensate are not co-rotating, under which condition, 
dynamical phases and phase transitions may exist [TtI ]. 
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Note added — During the preparation of the paper, we 
noticed the work of Reijnders and Duine 0] who studied 
some aspects of this system using a variational approach. 
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